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1. Introduction 



Given the study of quark-gluon plasma at RHIC and LHC, there is a continuing need for 
understanding what QCD actually predicts for the various physical observables of interest, in 
case the system does reach thermodynamical equilibrium. This is a highly non-trivial task, 
due to the strong interactions present in QCD even at temperatures above a few hundred 
MeV. Some of the observables can in principle be determined from first principles using lattice 
methods, but at present systematic errors, related particularly to light dynamical quarks, 
might be substantial. Therefore, there is a demand for alternative tools, whose results could 
be compared with those from lattice simulations. 

Because of asymptotic freedom, a natural candidate for such a tool is perturbation the- 
ory. Alas, it suffers from serious infrared (IR) problems at high temperatures ^H]. These 
problems are of a rather specific kind, however: they concern directly only static, long-range 
gluons. This has lead to the general framework of dimensionally reduced effective field the- 
ories [31 ■ The idea is that there are still a number of computations which can be carried 
out reliably in perturbation theory, while the parts which cannot, can be treated with lattice 
simulations using an effective theory which is much simpler than the original one, all the 
quarks having been integrated out. 

During the last few years, the program of dimensional reduction has been applied to a large 
variety of physical quantities: in particular, to various gluonic correlation lengths 1^- JH] as 
well as the equation of state ^J-^l]. Comparisons with 4d lattice results in pure SU(2) 
and SU(3) Yang- Mills theories are summarised in [3 |H1 El El EI > and they suggest that 
dimensional reduction is applicable as soon as the temperature is somewhat above Tc, where 
Tc is the critical temperature for the deconfinement phase transition. 

An even simpler pattern can perhaps be argued for jTJlElEl: it appears, inspecting the 
empirical results from the studies mentioned, that a truncated coupling constant series may 
already give a reasonable part of the full answer, if carried out to a sufficient depth such 
that all dynamical scales in the system have made their entrances. Therefore, either non- 
perturbative studies with the dimensionally reduced theory, or the further reduction suggested 
by this prescription, may provide a welcome alternative to direct 4d lattice simulations of 
QCD, particularly for observables where the systematic errors are not well under control. It 
should be mentioned that other proposals for analytic tools exist, as well (for recent reviews, 
see UHl). 

An important class of observables which have attracted somewhat less attention in these 
settings, are those built out of quark fields, rather than gluonic fields. In general, they are 
expected to be less IR sensitive and thus more perturbative than purely gluonic observables, 
so that analytic tools should be better applicable. At the same time, there is more urgency 
for analytic tools, given the severe difficulties in treating light dynamical fermions on the 
lattice. On top of these formal motivations, "mesonic" and "baryonic" observables may be 
more directly related to experimental signatures than gluonic ones. During the last few years. 
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significant progress has been made on one class of observables of this type, that is various 
quark number susceptibihties (see, e.g., The purpose of this paper is to address 

another class, closely related but more directly sensitive to IR physics: the long-distance 
spatial correlation lengths related to mesonic operators (i.e., gauge-invariant quark bilinears) 
constructed out of the light quark flavours ^2F. 

There is a large variety of independent mesonic observables of the type mentioned. Assum- 
ing, as we will in this paper, QCD to consist of Nf flavours of degenerate quarks, they can 
be divided, first of all, into singlets and non-singlets under the flavour group. Each class can 
furthermore be divided into scalar, pseudoscalar, vector, and axial vector objects. Physically 
perhaps the most interesting quark bilinears are those representing baryon number density 
and electric charge density, but since all the other ones can in principle be measured on the 
lattice equally well and this is the most immediate reference point for our results, we will try 
to be as general as possible. 

Previous analytic work on similar topics was initiated more than a decade ago |26j-|28j. 
following the first numerical studies [25]. (There was even earlier work for small temperatures, 
< T <C Tc [2nilSni-) Many of these studies had mostly qualitative motivations, however: 
corrections of order unity with respect to the ones computed were (knowingly) left out. They 
were also partly motivated by some qualitative features in the early lattice results (particularly 
that the vector channel, or p, appeared to be significantly "heavier" than the pseudoscalar 
channel, or tt, even at T > Tc), which have turned out to be finite lattice spacing artifacts. 
Significant further progress was made in j31j . and we show here, following that line, that the 
corrections neglected are numerically important, and even modify the parametric form of the 
final result, by removing any logarithms. Note that a somewhat related problem has been 
addressed in the context of the electroweak theory at finite temperatures |32j . although for 
operators which are bilinear in the Higgs field rather than in fermionic fields. 

The plan of this paper is the following. The basic notation used throughout is fixed 
in Sec. El Some well-known properties of mesonic correlators at very high temperatures are 
reviewed in Sec. El In Sec. 01 we discuss how the mesonic correlators can be incorporated into 
the general setting of describing high-temperature QCD by a dimensionally reduced effective 
field theory I31j . The quark part of the Lagrangian entering will be called three-dimensional 
(3d) "Non-Relativistic QCD", NRQCD3. The basic matching computation between the orig- 
inal QCD and NRQCD3 is carried out in Sec. El and the first non-trivial correction to the 
correlation length of flavour non-singlet mesons is then determined from NRQCD3 in Sec. 
In Sec. 13 we compare with lattice results. Flavour singlet observables, for which the pat- 
tern is qualitatively different from the non-singlets, since they may mix with purely gluonic 
operators, are treated in Sec. |H1 and we conclude in Sec. 
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2. Notation 



The theory we consider in this paper is QCD at a finite temperature T, with the Euchdean 
quark Lagrangian 

4 = V^(7;.I?^ + Af)V^ , (2.1) 

where = d^ — igA^, = A^^T'^, are Hermitean generators of SU(A'^c) normalised such 
that Tr T'^T'^ = \ 5"'^, and M is the mass matrix. For simplicity, we take M to be diagonal 
and degenerate, M = diag(m, m), and in most of what follows, m = 0. The number of 
flavours appearing in tp is denoted by Nf. 

Given tp, ip we may define, as usual, various bilinears. The Nf x Nf flavour basis is generated 

by 

^a_|^s^^n|^ i7- = l^^,^^^, n = l,...,Nf-l, (2.2) 

where F** = iNfxNf is an Nj x Nf identity matrix ("singlet"), and the traceless ("non- 
singlet") are assumed normalised such that 

Tr [F™F"] = i 5'"" . (2.3) 

We may then consider scalar, pseudoscalar, vector, and axial vector objects, 

S" = i)F°-ij , (2.4) 

P<^ = V575i^"V', (2.5) 

^ i^-i^F-tl^, (2.6) 

Al ^ V^T^Tsi^^. (2.7) 

The operators just defined may be interpreted either directly as physical condensates or 
currents, or as interpolating operators for certain particle states. In three-flavour QCD 
(A'^^ = 3) with a degenerate mass matrix, for instance, some possible assignments are shown 
in Tabled 

Given such operators, the correlators to be considered are of the form 

Cq[0", O^] = j^J^ dr j d^x e''^-'' (0"(r, x)O^(0, 0)) . (2.8) 

The corresponding configuration space correlator is 

C^[0-, O"] ^ J ^ e-'^-^ C7q[0^ O"] = fj^dr {0'^{t, x)O''(0, 0)) . (2.9) 

The expectation values are taken in a volume which is infinite in the spatial directions but 
finite in the temporal direction, with extent 1/T. Gauge fields are assumed to obey periodic, 
fermions anti-periodic boundary conditions around the temporal direction. 
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Operator Particle Behaviour in free theory Physical interpretation 





k/cj 




gn 


ao/(5 




ps 


•q' 


q^B3d{2po) 


pn 


TT,r] 


g253d(2po) 


V? 


to 


(g2+4p2)i?3d(2po) 


VI 




((72-4p2)i?3d(2po) 


T/n 


p,4> 


(g2 + 4p2)i?3d(2po) 


VI 




(g'-4pg)i?3d(2po) 


Ah 


fl 


(g2 + 4p2)53d(2po) 


Al 




(g2-4p2)B3d(2po) 


An 
^0 


ai 


(g2 + 4p2)B3d(2po) 


AI 




(g2-4p2)i?3d(2po) 



baryon density 
charge density 



Table 1: Different particle assignments and the nature of the singularity structure for the 
mesonic correlators considered (cf. Sec. The physical interpretation of Iq" applies for 
Nf = 3. The notation refers to polarizations transverse with respect to the spatial 

direction in which the correlation is measured; the longitudinal correlators all vanish in 
perturbation theory, due to current conservation in the chiral limit. 



The general structures of the correlators so defined are, because of rotational invariance, 
Cq[S",5''],Cq[yj^,^o'] = S'^'fiq'), (2.10) 

c^[vf,v^''] = ^'''[{^^,-'^)t{e) + ''-fKa')\. (2.11) 

and correspondingly with S ^ P,V ^ A. Here q= |q|. For conserved currents, l{q'^) = 0. 

We expect, in general, that in the full theory the analytic structures of the functions 
f {q^) ,t{q^) ,l{q^) around the origin are characterised by simple poles. The reason is that, as we 
will recall presently, these functions can after analytic continuation be viewed as propagators 
of physical states in a (2+l)-dimensional confining theory. The spectrum of the confining 
theory consists of bound states, represented by poles. Consequently, the Fourier transforms 
should show exponential falloff at large distances. The purpose of this paper is to determine 
the pole locations or, equivalently, the long-distance exponential falloffs of the corresponding 
configuration space correlation functions, to order 0{g^). The coefficients of the exponential 
falloffs are often referred to as screening masses. 
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3. Correlators at very high temperatures 



To start with, it may be useful to recall the behaviour of the correlation functions we are 
interested in at very high temperatures. In this limit, asymptotic freedom should guarantee 
that we can use perturbation theory, and the correlators can be computed with free fermions 
(Fig. ^a)). One finds, in dimensional regularisation, 

(3.1) 

where is the Dirac matrix appearing in O", = {1,75,7^,7/^75}, ^ = "JfiP^, and p„ 
denotes the fermionic Matsubara frequencies, p„ = 27rT(n + ^ ). The Dirac algebra and the 
integration are trivially carried out: the result contains either constants, which correspond 
to contact terms ~ 6{x) after going into configuration space according to Eq. 1)2. 9|) . or the 
function 

d^p 1 i 2pn — iq 



BM{2pn) - J (2^)3 + p2][p2 + (p + q)2] " 8nq 2p„ + " ^^'^^ 

The qualitative functional forms of the terms including this function for the different are 
shown in Table Q for p„ = Po- 

The expression in Eq. (|3.1j) contains a sum over all pn = 27rT(n + ^ ). Since the non- 
trivial structure appears at the point 2p„ (cf. Eq. (|3.2|) '). the correlator is dominated at large 
distances (small q) by the smallest Matsubara frequencies, itpo = ±7rT. These frequency 
modes live in three spatial dimensions. If we now imagine viewing the three-dimensional 
theory rather as a (2-|-l)-dimensional one, and call the direction in which the correlation is 
measured, the time t, then the screening mass corresponds to the energy of a state consisting 
of two on-shell massive quarks, each of "mass" pq. Put in this language, our problem will be 
to compute the first correction to the energy of the two-particle state. 

We may also note that the nature of the singularity at q = ±2ipQ in Eq. 1)3. 2|) is a branch 
cut. In a (2-|-l)-dimensional world, this just corresponds to an energy threshold for producing 
two free on-shell particles of "mass" po- We may expect this structure to get qualitatively 
modified by the interactions, which place the quarks off-shell and bind them together, as just 
described: the threshold should then convert to a pole. The computation of the functional 
form of the singularity is beyond the scope of the present paper, however: we simply study 
how the location of the threshold at 2pQ gets modified. 

Note from Table ^ that for the charge-charge correlators (V^'^,^q), there is a prefactor 
(g^ -|- 4pq) which vanishes at the point of the branch singularity in i?3cj(2po)- This means 
that a contribution of the singularity at g = ±2ipo is suppressed, and the correlator tends 
to decay even faster: it is easy to see that in configuration space, it is power-suppressed in 
1/pot. The issue is even clearer if, for a moment, one imagines going to a world with one 
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(a) 



(b) 



(c) 



Figure 1: (a): The leading order correlator, (h), (c): Graphs representing non-trivial gluon 
condensates in the operator product expansion for the mesonic correlation function. 



space dimension only: then 
Sid(2po) 



dp 



1 



PO g2 + 4p2 ' 



(3.3) 



2t: [pl+p^][pl + {p + qf] 

and the singularity (a real pole in this case) would be completely cancelled by the prefactor, 
leaving only (short-distance) contact terms^. We will return to this issue below. 

To conclude this Section, we briefly inspect the correlators from yet another angle. Let us, 
for a moment, consider the correlator not around the pole at g = ib2ipo; but for a very large 
(real) q. Then the correlator can be computed in the operator product expansion: taking a 
background of constant (zero Matsubara mode) gauge fields (Figs, ^b), (c)) one obtains, for 
instance, 



Nr 



2po - iq 



Apo + iqln 
47r V 2po + tq 



+ 



1 



Trpo (g2 + 4p§)2 



g^TrAl + ... 



(3.4) 



Extrapolating then back towards the threshold region, say q ~ 2ipQ + 0{g'^T), one sees how 
the corrections related to the gauge fields are becoming increasingly important. In order to 
systematically evaluate them we have to turn, however, to other methods. 



4. Basic structure of NRQCDg 

To go beyond the accuracy of Fig. ^a), we need to compute graphs of the qualitative types 
in Fig. 121 Just the graphs shown are not enough, however: one also needs to account for 
various kinds of iterations of these basic topologies. The reason is that near the threshold, 
quarks can be almost on-shell, |1/^| ~ 0{l/g^T), and thus compensate for the appearance 
of 0{g'^T) in the vertices. Therefore, there is a large set of graphs which has to be identified 
and resummed, in order to obtain the first correction to the threshold location. 

As in many other contexts, a convenient organising principle for carrying out such resum- 
mations is offered by effective field theory methods. As already explained, in the present 

^This remains true even after inclusion of the other pn's. 
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(a) (b) (c) 

Figure 2: Schematic representations of various classes of beyond-the-leading order corrections: 
(a) quark self-energy correction; (b) quark-antiquark interaction via gluon exchange; (c) for 
flavour singlets, the correlation can also he mediated by purely gluonic states. 

case we can view the correlation lengths as (2+l)-dimensional bound states of heavy parti- 
cles of "mass" poi rnuch larger than the infrared scales ~ gT,g'^T of gauge field dynamics 
at high temperatures P]. The masses of such heavy bound states can be determined, like 
those of quarkonia at T = 0, through an effective theory called the non-relativistic QCD 
(NRQCD) 33 . Our case does lead to some significant differences with respect to this frame- 
work; for instance, that the "masses", po, are not really masses in the conventional sense, 
but conserve chiral invariance and are directly scale invariant physical quantities. Despite 
the differences, we will refer to the general framework to be used as NRQCD3. The philos- 
ophy of applying NRQCD techniques to the present problem was introduced by Huang and 
Lissia |31| . We start by showing how the non-relativistic structure emerges on the tree- level. 



4.1. Tree- level 

As seen in Eq. ()3.1() . at leading order the correlators considered get independent contributions 
from the various quark Matsubara modes, with the lowest ones being dominant at large 
distances. Let us therefore again concentrate on pQ = ttT. The quark Lagrangian for this 
Matsubara mode is 

C% = %1} i-foPo - igioAo + ikDk + 73D3 ifj , (4.1) 

where A; = 1, 2. 

In Eq. (|4.1j) , ip interacts with gluonic zero Matsubara modes only. The Lagrangian describ- 
ing their dynamics, obtained after integrating out all quarks and the non-zero Matsubara 
modes of gluons, is the dimensionally reduced Lagrangian, of the form [2j 



£^ = i Tr 4 + TV [A, Ao]' + m|Tr + X^^\Tr A^f + A^^Tr + ... . (4.2) 



Here ^ = 1, 3, Fij = {i/gE)[Di, Dj], and Di = di — igE^i, where = A^^T^ and the fields 
have now the dimension GeV^/^, after a trivial rescaling with T^/^. The parameters, to the 
orders that we need them, are 



m|=(^ + ^)5^^^ 9l = g^T, (4.3) 
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while X^^^ = 0{g^T) can be ignored. The parameters have been determined up to next-to- 
leading order P^^nEHl- According to Eqs. ()4.2() . ()4.3() . infrared gauge field dynamics is 
sensitive to two scales: the perturbative scale ttt-e ~ gT associated with the colour-electric 
component Aq, and the (2-|-l)-dimensional confinement scale g"^ ~ g'^T associated with the 
colour-magnetic components Ai. 

Let us next set up some notation. In the following, we use rotational invariance to choose 
the vector q in Eq. (|2.8I) to point in the X3-direction. According to Eq. 1)2. 9|1 we are thus 
effectively considering correlators averaged over the (xi, a;2)-plane, and depending on X3. 
We will often denote the Euclidean xs-coordinate by t, and the vector (xi,X2) by x_l. The 
corresponding Fourier modes are ps, p_L. The indices fc, I are assumed to have the values 1, 2, 
and label the transverse directions. Given this choice of q we can consider, without loss of 
generality, the correlator 



Ct[0\0' 



dga 
2tt 



l/T 



dr /dV(O"(r,x^,t)O''(0,0x,0)) , 



(4.4) 

rather than Cx[0",0^] in Eq. (tToll . 

Once this choice has been made, the quark Lagrangian in Eq. (|4.1|) can be written in a 
suggestive form. Let us make a basis transformation from the standard representation 



70 






where are the Pauli matrices, to a new basis, 7" 




75 



[l^standardf/-!^ with 



(4.5) 
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1 







1 





1 
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1 





-1 
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In this basis, 

, / (po + ^P3)l2x2 -ekiPkcri 
where eti is antisymmetric and 612 = 1- Thus, denoting 



(4.6) 



(4.7) 



X 



where 4' s-re two-component spinors, the Lagrangian in Eq. (|4.1j) becomes 



E = ^X 



Po - 9A0 + D3 



X + i(, 



Po - 9A0 - D3 + (/)hkiDkaix - X^ekiDk 



(4.9) 
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Leaving out, for brevity, the flavour structures, the biUnear operators of Eqs. H2.4p ~ H2.7|l 
become 



S 
P 
Vo 
Vk 

^3 
^0 

Ak 



-^(x"^o-3X + (A^CTS^) • 



'lp{T, X) 



7oV'(t, -x), reads x(t,x) 



(4.10) 
(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 

, -x). 



Note that the parity transformation, 
(/)(r, x) x{tj ~x) in this basis. 

In Sec. IHl we saw that in the correlators we are interested in, the quarks are almost "on- 
shell". Let us therefore consider Eq. (|4.9)) for on-sheh configurations, in momentum space. 
We fix and consider the functional dependence on p^. In the free case, the on-shell point 
is at Po + p'^ = 0, i.e., p^ = ±i[po + pi/(2po) + •••]• The positive and negative "energy" 
solutions could be called the particles and the anti-particles. 

Because the quarks interact (at the tree-level) with bosonic Matsubara zero modes only, 
we expect that off-shellness is related to the momentum scales of the latter: {p^ ± ipol <gT. 
Assuming this to be the case, the idea is to construct an effective theory describing the soft 
scale dynamics. For instance, consider the pole around ~ ipQ. Then x represents the 
light mode, while is heavy. Being so far at the classical level, we can solve for cj) from 
the equations of motion, and substitute the solution back to the Lagrangian, resulting in an 
action for x alone (cf., e.g., the pedagogic presentation in for the 4d zero temperature 
case). The same treatment naturally applies to (j) around its pole. Expanding finally in 
powers of 1/pO) we obtain a "diagonalised" on-shell effective Lagrangian for two independent 
light modes, with a non-relativistic structure: 



iX 



+ i4>^ 



Po 



Po 



gAo + Ds 
gAo - Ds 



2po 



X 



9 



+ 



(4.18) 



The mesonic states are most usefully still represented by the operators in Eqs. 1)4. 1U() " (|4. 17(1 . 
For future reference, let us note that the free propagators following from Eq. (|4.18|) are 



{Xu{p)x*v{(l)) 

{MpWviQ)) 



5uv{27rf6^^\p-q) 
6uv{27rf6^^^{p-q) 



M + ip3+ pi/(2po) ' 
—i 

M-ip3 + pi/(2po) ' 



(4.19) 
(4.20) 
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where M = po, and u,v contain the spinor, flavour and colour indices. We may also remark 
that going into configuration space, the propagators become 

{Xuix)xliy)) = -iSu.eixs-y^):^ -exp(-M|x3-^/3|- ^°^f^~^V' )' (4-21) 

(M^my)) = -i6uve{y3-X3)—^ -exp(-M|x3-y3|- ^°^''^~^V' )- (4-22) 

27r|j;3-y3| V 2|j;3 - ^3] / 

This implies forward time propagation for x and backward for (p. For pQ ^ 00, 

\2 



PO 



V 2\xt — / 



2vr|x3-?/3| V 2\x3-y3 

indicating that the quarks become static. 

To conclude this Section, it may be noted that in Eqs. H4.1U() - ()4.17() the charges {Vo,Aq) 
and the third components of the currents (V3, ^3) are of the type ~ x^X + 4'^4'- This implies 
that the corresponding correlators decay very rapidly: according to Eqs. (|4.2H) . (|4.22|) . they 
in fact vanish identically at non-zero distances, being of the form ~ 9{t)6{—t). For the third 
components of the currents this is due to current conservation (the correlator is a constant, 
and vanishes in infinite volume, up to possible contact terms), while for the charges, this is 
related to the discussion in Sec. |21 where we observed the prefactor [q^ + Ap^) in the free 
correlators. In the on-shell limit, this causes the correlators to vanish; in reality, they are 
power-suppressed in l/p^t. In the following, we will concentrate on the correlators which 
decay more slowly, that is, which are represented by operators of the type ~ x^(/> + (j)^x- 



4.2. Power-counting and loop corrections 

The construction of the effective theory in Eq. (|4.18|) was so far on the tree-level. In order to 
be sure that the IR structure of the effective theory is really that of finite temperature QCD, 
and in order to determine the radiative corrections we are ultimately interested in, we need 
to consider loop corrections. In particular, we should find out the order to which one needs to 
go in the matching computation between the original finite temperature QCD and NRQCD3, 
as well as the resolution with which the dynamics inside NRQCD3 needs to be treated, given 
that we are interested in the 0{g^) correction to the spatial mesonic correlation lengths. For 
these purposes, it is useful to set up some parametric power-counting rules. 

Since the general upshot will be rather obvious, let us state it immediately: in order to 
determine the 0{g'^) correction to the mesonic correlation lengths, the additive parameter po 
in Eq. (|4.18|) becomes a matching coefficient (to be denoted by M, as in Eqs. H4.19() - (|4.22() ) 
which has to be determined to 1-loop order. On the other hand, no other parameters need 
to be matched at this level, and some interactions can even be dropped out: the final form of 
the theory is displayed in Eq. (|4.25|) below. This matching takes care of the "hard part" in 
loop integrals of the type in Fig.Efa), while the dynamics within the effective theory accounts 
for diagrams of the type in Fig. Efb), where the gauge fields are always soft. 
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To account for the full correlator defined in Eq. ()2.8|) . one should also carry out operator 
matching to order 0{g^). We will not need to consider it here, however, since we will not 
attempt to determine the overall magnitude of the correlator, only the coefficient of its 
exponential falloff. 

Let us now turn to the power-counting rules. As mentioned, the parametric scales in the 
problem are ~ ttT, gT, g'^T. Let us write m = TrT,v = g, and denote by Ap^ the off-shellness 
of the "energy", Ap^ = it ipQ. Following the terminology of NRQCD, one can define four 
different regions of phase space: 



hard (h) 
soft (s) 
potential (p) 
ultrasoft (us) 



Ap3 ~ |p_l| ~ m , 
Aps ~ |p_l| ~ mv , 
Ap3 ~ mv^ , I p_L I ~ mv , 

Aps ~ |p_l| ~ mv^ . (4.24) 



By definition, the NRQCD3 theory in Eqs. (|4.18jl . (|4.2jl is chosen to contain at most soft 
quarks and gluons as dynamical degrees of freedom. In fact, close to on-shell, the quarks of 
NRQCD3 are potential, since Ap^ ~ |pj_|/m ~ mv"^. 

In order to estimate the importance of the various momentum regions in Eq. (|4.24|) . it 
is helpful to assign power-counting rules also to the corresponding field components. This 
can be done by going to configuration space and requiring that the action, a pure number, 
is parametrically of order unity. For soft and potential quarks, Xs^Xp; it follows from S ~ 
/ dt d^x^ x^^f^tX ~ 1 that x ~ l/|x±| ~ mv. For soft gluons As, S ~ Jdtd'^:x.± Asd^Ag ~ 1 
implies that Ag ~ (t/x^)^/^ ~ m^/^u^/^. For ultrasoft gluons vlus; correspondingly, A^s ~ 
(t/xj_)^/^ ~ m^/'^v. In addition, when operating on on-shell quarks, the time derivative can 
be estimated as 9f ~ 1/i ~ mv"^ . 

Given these rules and that g(| ~ mu^, we note that ^e^s ~ mt>^/^, ^e^us ~ mv"^. The 
important conclusion follows that soft and ultrasoft gauge fields are parametrically of higher 
order than the transverse momentum ~ mv. In particular, given that we are interested in 
quarks which are off-shell by at most 0{mv'^) ~ 0{g^T), we should note that this magnitude 
is already saturated by V^/m: no transverse covariant derivatives are needed in the quark 
action! On the contrary, the gauge fields are significant in relation to dt. 

To summarise, if we want to keep the Lagrangian up to and including the order 0{m^v^), 
such that we know the on-shell quark self-energy up to order 0{mv'^), it is enough to re- 
place Eq. (|4.18|) with the action 

= ix^ (^M - 5E^o + - ^) X + - 5E^o - A - ^) , (4.25) 

where Dt = = — ig-^A^, and the gauge fields have now the same 3d normalisation as 
in Eq. ()4.2() . Furthermore, the only parameter requiring a matching computation is M. 

To conclude this Section, let us briefiy view the dynamics of the remaining degrees of 
freedom. Since the quarks are very heavy, of mass po, we may compute the static potential 
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in (2+1) dimensions. Its structure is, 

V{r) ^ gl\nr + g^r + 0{g%r^) 



(4.26) 



If inserted into a Schrodinger equation, we get 

^l^~^W~5|lnr, (4.27) 



and therefore |p^| ~ d/dr ~ 1/r ~ y g^po ~ gT. We see that the expansion parameter 
related to the static potential is 0{g\r) ~ 0{g). Consequently, for the correction of order 
E ~ |p±P/po ~ g^T, it is sufficient to compute the static potential to leading order. The 
leading order contribution itself may, however, already be infrared divergent, in the sense 
that the scale appearing inside the logarithm in Eq. (|4.26j) remains to be determined. 



5. Matching from QCD to NRQCD3 

As explained above, the parameter we need to determine through a matching computation 
is M, as defined by Eq. ()4.25() . This computation was not carried out in j^J, although many 
other matchings related to NRQCD3 were. 

In order to avoid having to deal with any wave function normalisation factors, we will 
carry out the matching by computing the quark "pole mass" , or finite temperature Euclidean 
dispersion relation, both in the original finite temperature QCD, as well as with the theory 
of Eq. 1)4. 25() . Both computations produce a gauge-invariant result, and requiring the out- 
comes to be equivalent, leads to an expression for M in terms of the parameters of finite 
temperature QCD. 

In order to monitor explicitly the gauge parameter dependence and the IR sensitivity of 
the computations carried out, we take the gluon propagator to be of the form 

{A1{p)Al{q)) = 6'^\27,f-^^5^'-^^\p + q) 

That is, we introduce a gauge parameter ^ as well as an IR-cutoff A. As we will see, the pole 
masses, computed in both ways, are independent of ^, and allow to take the limit A — > 0. 
The same is then also true for their difference, determining the matching coefficient we are 
interested in. 

An important and somewhat non-trivial point is the way the computation is to be carried 
out on the side of NRQCD3. The "problem" is that we want to use dimensional regularisation, 
yet obey the power-counting conventions as discussed above. If one, however, takes the free 
propagators as determined by Eq. (|4.25j) . shown in Eqs. (|4.19|) . (|4.2n|) . then (irrespective of 
whether or not one keeps M there explicitly or shifts it away by modifying the integration 
contour for ^3) momenta of order |p_l| ~ po do contribute. This, however, we do not want. 



p2 



+ A2 



+ 



P^,Pu 4 

p2 p2 _|_ e_^2 



(5.1) 
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because our power-counting rules assumed the dynamical scales inside the effective theory to 
be at most soft, |p_l| <mv, while po ^ m. The way to avoid the problem with dimensional 
regularisation has been clarified in 37-;,: one should carry out the matching in the soft region 
of the phase space, ps ~ |p^| ~ mv, where the transverse part p5_/po can be treated as a 
perturbation. That is, we match the operators of NRQCD3 order by order in 1/po- 

The situation would be different if we were to treat NRQCD3 also on the lattice. In 
order to match from dimensional regularisation to lattice regularisation, say, one would then 
really have to compute using the same forms of propagators in both schemes. This leads to 
additional logarithmic ultraviolet divergences in the parameter M, needed in order to cancel 
the corresponding divergences appearing in the dynamics within NRQCD3. Carrying out the 
computation this way is not much more complicated than what we do here, but since the 
results are not needed for our present purposes, we do not display them. 

On the side of QCD, the inverse propagator evaluated at the position of the tree-level pole, 
= 0, reads 

where (..•)/ denotes fermionic Matsubara modes, (...);, bosonic, and Cp = {N^ — l)/2Nc. 
The result comes from the transverse part of the gluon propagator; the longitudinal gauge 
parameter dependent part does not contribute at the pole. The task is now to find the zeros 
of to relative order 0{g'^). 

In order to streamline the task, we may set = 0, and multiply from the left 

with 7o. The result obtains a simple form if we choose the basis of 7-matrices introduced 
in Sec. 14.11 Indeed, off-diagonal terms vanish after integration over q, and the result becomes 
block-diagonal. Moreover, the upper two components are degenerate, as well as the lower 
two components, and the latter follow from the former simply by (j»3 — > —^3,^3 — > —^3). 
Therefore it is enough to consider 



(5.2) 





7oS(p)] = i{po + ip3 - g'^Cp^-, 2 I \2\ \+Po - '^P3 -Qo + iqs (Aq) (5.3) 

-po + ip3 + qo- m (M) 

-2po - 2ip3 + 2qo + 2iq3 (Ak)] } , 



where the gauge field components leading to the various contributions have also been indi- 
cated, and we have set the spatial dimensionality to d = 3, anticipating the fact that the result 
is UV-finite. The 0{g^) term is to be evaluated at the tree-level pole position po = —ip^. 

It can be observed from Eq. (|5.3|) that both ^0 and A2, give possibly logarithmically IR- 
divergent integrals (if A ^ 0), from the terms with pQ — ip^ in the numerator, but the sum 
is finite since they come with opposite signs. For A^, on the other hand, the contribution 
involving comes in the combination po + ^P3; and vanishes altogether at the pole. The 

correction as a whole is therefore IR finite, and can be evaluated with A — > 0, whereby the 
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integral becomes 
J: 2(go + ig3) 

Consequently, 



1 ^ 2p ■ q 



Po=-Ws 



Po%{p - (l)){(p)b 
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and the dispersion relation, or pole position, becomes 



P3 



rp2 

Po + 9^Cf — 
8po 



(5.4) 



(5.5) 



(5.6) 



On the side of NRQCD3 we proceed as explained above, following [^Tj. The gauge field 
propagator is treated as in Eq. ()5.1() at the current stage (during matching, any possible IR 
cutoffs should be the same in the two theories). The inverse propagator now becomes 

d^-^'q 1 



(27r)3~2f M + ip3- iqs 



+ 



1 

1 



(^0; 
(^3) 



(5.7) 



g2 + A2 

The potential IR divergences from the two gauge field components are the same as in Eq. (|5.3j) . 
but they cancel, as they did there, such that the 1-loop correction is finite, and in fact vanishes. 
Thus, the dispersion relation is simply 



Ps^iM . 

Matching now Eqs. (|5.6|1 and (|5.8|1 . we obtain 

M = po + g^T^ + 0{g^T). 
Sir 



(5.8) 



(5.9) 



6. Perturbative solution within NRQCD3 

After the construction of the effective theory, which accounts for hard gluons in Fig. ISfa), 
we still have to account for the soft part of Fig. Ufa) as well as for all of Fig. I^fb), where the 
gluons can only be soft (i.e., Matsubara zero modes). This is to be done by computing the 
analogues of the correlators in Eq. 1)4.4(1 . using the effective theory in Eqs. I|4.2() . (|4.25() . and 
the operators in Eqs. (|4.1U() - (|4.17() . Note that because of the appearance of Debye screening 
(represented by the parameter in Eq. ()4.2|) '). the gauge field components Aq^A^, behave 
now differently: the propagators are 

(A^oipXiQ)) = 6^\27Tf-'^S^'-'^\p + q)^^^ , (6.1) 

{Anp)A)iq)) = 5"^(27r)3-2^5(3-2e)(p + ^) 



PiPj 



p^ / p 



+ A2 



+ 



pipj ^ 

p2 r)2 + f A2 



(6.2) 
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We have again introduced a fictitious infrared regulator A, as well as a general gauge param- 
eter ^, just in order to check that the results are independent of them. 
The correlators of Eq. (|4.4j) take the form 

a[0^0''] ~ y'dV(o"(x±,t)O''(0x,0)) , (6.3) 

where we have suppressed the dependence on the r-coordinate which disappeared as we 
moved to Matsubara frequency space. As already mentioned, we are ignoring corrections of 
relative order 0{g'^) in the overall normalisation of Eq. (|6.3|) . but not in the coefficient of the 
exponential falloff. 

Even though the operators in Eqs. 1)4. lU() - ()4. 17(1 contain two terms, ~ x^4^ + </'^X) only one 
product contributes in Eq. (|6.3j) : this is the one where both quarks are forward propagating 
(cf. Eqs. (lOTT) . (11221)), ~ ((/>t(t)x(t)xHO)(/>(0)). In other words, we need to determine the 
Green's function for the composite object / d^x^ 0* (x^, t)ru„Xt;(x±, t), where Tuv is a 2 x 2 
matrix in spinor space, an Nf x Nf matrix in flavour space, and an x Nc matrix in 
colour space. Since the spinor and flavour indices play no role in the following, they will be 
suppressed. 

As we have argued in Sec. 14.21 it should be rather obvious how to proceed now: simply 
integrate out the gauge fields, to produce a static potential for the (j)*x pair, and solve 
the corresponding Schrodinger equation. The reason this works is that, as we will show a 
posteriori, the non-perturbative ultrasoft gauge fields do not contribute at the present order. 
The problem could also be formalised for instance by going through the PNRQCD effective 
theory from which the soft gluonic scales have been systematically integrated out. 

For our modest purposes here, involving only the Coulomb-type potential, going through the 
full formalism of PNRQCD seems somewhat too elaborate, however. What we rather wish 
to do is to recall on a concrete but somewhat heuristic level how the Schrodinger equation 
can be seen to emerge, just in order to get all the signs right in our final expression. 

The Green's function we need to consider is Eq. ()6.3|1 . We would like to find the equation 
of motion obeyed by this Green's function, which we expect, at least at large t, to be of the 
form {dt — H)G{t) = C6{t), where H is a differential operator and C is a constant. Then the 
exponential decay is determined by the lowest eigenvalue of H. 

In order to find H, we now define a point-splitting (which is essentially gauge- invariant 
because no spatial components of gauge fields appear in Eq. (|4.25jl ). by introducing a vector 
r and, with the intention of taking r ^ afterwards, consider the modified correlator 

C{r,t) = y'd2-2^R(r(R + r/2,t)x(R-r/2,t) x*(0, 0)0(0, 0)) . (6.4) 

We denote the tree and 1-loop contributions by C^'^\r,t),C^^\r,t). It is then easy to see 
that 

dt + 2M- —Vl] C W (r, t) oc 5(t)(5(2-2^) (r) . (6.5) 
L Po J 
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The result for the 1-loop contribution (a sum of the first two graphs in Fig. ^ can, on the 
other hand, be written in the form 



dt + 2M - — V 
Po 



1 Vr- ml 



CW (r, t) = -glCpf^ ^ 7^' ^' ^ ' ^' ' (6.6) 

\^Po Po Po Po / 

where the kernel /C is dimensionless. Since the kernel is multiplied by g^, we may assume 
that it can be expanded with respect to its first two arguments, given that non-trivial terms 
in this expansion should be subleading in g^/po with respect to the terms dt, — V^/po already 
appearing in Eq. 1)6. 5|) . The zeroth order term turns out to involve a momentum space delta- 
function 5(^3), which restricts the remaining loop integral to be (2 — 2e)-dimensional, and 
produces just the 1-loop static potential. 



V{r) 



V Po Po 

d2-2^Q 



(2vr) 



2-2e 



1 



A2 



1 



1 



+ mi 



1 + e 



Writing now C{r,t) = C^°\r,t) + C^'^\r,t), we can combine Eqs. dHSl), to 

dt + 2M- —Vl + y (r)l C7(r, t) oc 5{t)S^^^^'^ (r) , 
Po J 



(6.7) 



(6.8) 



where the point-splitting can be removed, r ^ 0, after the solution, and subleading terms 
have been dropped. Therefore, the exponential falloff is determined by the eigenvalues of the 
Hamiltonian read from Eq. (|6.8j) . 

It is worthwhile to note, in passing, the signs appearing in the two terms in Eq. ()6.7() (the 
first from A3, the second from Ao). The reason for the differences are the ways in which 
A^jAo appear in Eq. 1)4. 25|) : ztig'E^a versus — (^e^o- The fact that the r- independent terms 
from A^^Ao come with opposite signs, was already met in Eqs. (|5.3() . 1)5. 7p . 

Employing the basic integrals 



d- 



2-2e, 



1 



-2e 



(2vr)2-2e q2 + ^2 
2^ ^jq-r 



- + 21n^ + 0(6) 
47r V e A 



(6.9) 

where p, is the MS scale parameter, fx = fl{e^ /Airy/'^ , and Ko is a modified Bessel function 
as well as the asymptotic behaviour of Ko, 



iKoiXr) , 



Ko{x) = -In--^E + O{x) , 



(6.11) 



it is observed that Eq. (|6.7|1 is both ultraviolet finite (lim£_>o exists) and infrared finite (lim^^^g 
exists). Its final form is 



Vir) 



2 Cp 



niEr 

In— \-jE- KoimEVj 



(6.12) 
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The potential determines the coefficient of the exponential fall-off, ^ 



-1 



2M 



^ + Vir) 
Po 



*0 = "IfuU^C 



mfuii, through 

(6.13) 



where ^'o is the ground state wave function. Therefore, mfuii is indeed insensitive to the glu- 
onic ultrasoft scale at this order: there are no ambiguities inside the logarithm in Eq. H4.26() . 

In order to determine nif^w, we have to solve Eq. (|6.13j) . We are not aware of an analytic 
solution for this bound state problem. In order to find the solution numerically, we rescale 



r 



mull 



2M 



n2 P 



whereby Eq. H6.13() becomes 



dr^ f df 



+ p( In - +7B 



where 



P 



PoglCp _ 3(iV2 - 1) 



27rm| 



2A^,(2A^, + Nf) 



(r)- 




^-0 = 


2/3 


,Nj 


= 


1/2 


,Nf 


= 2 


4/9 




= 3 



(6.14) 



(6.15) 



(6.16) 



A numerical solution of Eq. (|6.15() is facilitated by determining the analytic behaviour of 
^^{r) around the origin. Assuming ^'o(O) to be finite, one finds 



1 



1 



^'o(f) pa ^-0(0) l + -pr'[\n- + ^E-l--E, 



Integrating towards large f and requiring square integrability, gives then 



£"0 = 0.16368014, 


P = 


2/3, 




= 0) 


= 0.38237416, 


P = 


1/2, 


{Nf 


= 2) 


= 0.46939139, 


P = 


4/9 . 


{Nf 


= 3) 



(6.17) 



(6.18) 
(6.19) 
(6.20) 



Therefore, for Nc = 3, Nf = 
(|6.18j) together imply that 

mfuu = 27rr 



(corresponding to quenched QCD), Eqs. (|5.9|) . (|6.14jl and 
! 27rr + 0.14083730 5^ T . (6.21) 



2 / I 
^^2^12 



Er 



It is worthwhile to note that this expression contains no logarithms. Here we differ from |26j . 
for example, where the contributions from were left out. As can be seen in Eqs. (|6.7j) . 
()6.12() . these contributions are essential. 

At T ~ 2Tc for Nf = 0, it is estimated jSj that the effective gauge coupling is g'^/T ^ 2.7. 
Thus the second term in Eq. (|6.21j) represents a ~ 5% correction to the leading order result, 
with a positive sign. Since the spinor structure of the operators in Eqs. (|4.1()j) - (|4.17j) played no 
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significant role in the analysis above (it only leads to an overall factor, Tr 1 or Tr o"!), the same 
result holds for all the flavour non-singlet operators of the type (j)^Xi that is, S"", P", Vj^, A^. 
As we have mentioned, the correlators related to Vq'jAq are suppressed, while those related 
to V^jA^ are constants because of current conservation, and vanish in infinite volume, such 
that no meaningful long-distance screening masses can be associated with them. 

7. Comparison with lattice 

Let us now consider lattice measurements of the correlation lengths discussed, at T > Tc. This 
is a classic problem in lattice QCD, initiated already a long time ago 25 . Most measurements 
concentrate on the non-singlets P", (or "vr, p"), because the singlets require an additional 
"disconnected" contraction (corresponding to the quark lines in Fig. [2t^c)), whose inclusion 
is non-trivial. For recent original work on the topic, we refer to j4()j-j42j. and for a recent 
review, to (16) . 

As always with lattice simulations, practical measurements are faced with finite volume 
and discretisation artifacts. The measurements we are concerned with face also further prob- 
lems, since the operators are fermionic: we would like to approach the chiral limit and 
simultaneously respect the flavour symmetries, which turn out to be hard requirements to 
satisfy. Simulations could be carried out with, say, staggered, Wilson, domain wall, or over- 
lap quarks, but in all cases but the last, which is very expensive in practice, it is oftentimes 
difficult to estimate how serious the systematic errors are. Finally, most lattice simulations 
are "quenched", ignoring the fermion determinant, which we however do not need to consider 
a restriction, since we can easily set Nf = in our predictions. 

To summarise the current status, from quenched measurements with Wilson quarks at 
Nj- = 8 |1D]: already at T ~ 1.5Tc, the tt is just slightly below p. Both have masses consistent 
with or very slightly below the free theory prediction, that is 2ttT in the continuum. Very 
close to Tc the situation has to be different, of course, and the pion mass should even vanish 
at Tc while approaching the chiral limit with dynamical quarks, for Nf = 2. 

This pattern is quite different from what it used to be: in the early days, staggered quarks 
gave an anomalously low screening mass for the pion even at T > Tc, both in the quenched as 
well as in the unquenched case. It has been found out, however, that discretisation effects are 
substantial for staggered quarks: in the quenched theory, the light masses increase in going 
from A'^^ = 4 to At- = 12, at T = (1.5...2.0)Tc, by up to ~ 40% (when masses are expressed 
in units of T) j4H \2'A\ . On the other hand, the effects of quenching seem to be moderate, 
although only tested far from the continuum limit (for A^^ = 4, staggered quarks) j42j . 

To summarise, within the current resolution the mesonic screening masses agree with the 
leading order prediction, within statistical and certainly within systematic errors. It will 
require a significantly higher resolution to show unambiguously whether the 5% correction 
we have found is really reproduced by the lattice data in the chiral continuum limit. In any 
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case, the classic problem, that the vr and p seemed not to be degenerate sXT > Tc in contrast 
to the perturbative prediction, has disappeared while approaching the continuum limit. 



8. Flavour singlets 

Flavour singlet quark bilinears behave in a different way than flavour non-singlets, because 
their correlation functions contain a "disconnected" contribution, of the type in Fig. Efc). 
This means that the operators couple, in general, to bosonic glueball states. They are also 
more difficult to measure on the lattice, just because of the disconnected contraction. Here 
we discuss some basic features concerning the flavour singlet correlators. 

The relations of the flavour singlets to gluonic operators can be determined by coupling the 
flavour singlets to sources in the Lagrangian, and then performing the Grassmann integral 
over the quarks, in a given gauge fleld background. At high temperatures, this procedure 
is infrared safe, and the results can be expanded in m^/(7rT)^. To leading order in the 
expansion, a 1-loop computation of 0-gluon, 1-gluon, 2-gluon and 3-gluon vertices leads to 

^ Nj^-^g^TiAl, (8.1) 

_L^^o^ ^ -^g^TvAl, (8.3) 

^ 0. (8.4) 

These relations imply, for instance, that the correlation length of S'^ = i/jip can be determined 
from Tr^Q, and the correlation length of the baryon density Vq = "ip^Qip from Tr^Q [9} I44j . 
Furthermore, from 

;(^al3-y5Fal3^yS = d^K^ , (8.5) 



327r2 °'f^'y^ '^1^ 1^ 
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it follows by integrating over J dr d2x_|_ that the correlation length of = tpj^ip is determined 
by that of oc g'^eij^AQF-j, where the form indicated applies in the static limit. All of these 
gluonic correlation lengths have been measured, for various A'^^, in ref. jH]. It may also be 
noted from Eqs. (|8.1() . ()8.2() that some of the couplings vanish in the chiral limit. 

For the axial charge Af^ = V'7o75'0) relation to gluonic objects is very non-trivial because 
of the axial anomaly. We recall here only that if one couples the axial charge to a constant 
chemical potential, and then integrates over spacetime, one obtains [IH] 

-^V'ToTs^ ^ Kq ■ (8.7) 
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As is well-known this relation does not unambiguously apply for local operators, however (the 
right-hand-side, the Chern-Simons number density, is not even gauge-invariant). Nevertheless 
it shows explicitly the quantum numbers of the gluonic objects that may couple to. 

Let us finally discuss the spatial components of the vector and axial currents, and also 
touch the issue of Ua(1) axial symmetry non-restoration (for a review see, e.g., |1Z|)- To 
start with, we may recall that 'tpj/^'J^ip and 'ipj^ip are related through the anomaly equation 

d^[^P7f,75^] = 2m V'TsV' + N f^e^p^sF^pF^s ■ (8-8) 

Indeed, this relation remains true even at finite temperatures |45j, in spite of the fact that 
the ways in which 9^ ["07^ 75 V'] and ip^^ip separately couple to gluonic operators, change 
Integrating now (the chiral m ^ limit of) Eq. 1)8.8(1 over /drd^x^ and noticing that the 
partial derivatives ^3 appearing on both sides do not modify the qualitative structure and 
can thus be dropped out, it is observed that the correlators in the X3-direction of ^3 = 
''/'7375V' are determined by the operator oc g^eij^A^F^j. Thus the correlator of ^3 falls 
off exponentially, with a correlation length determined |481 E] to be ^"^[^3] = ^-[^3] = 
m[i^3] niE + g^Nc{ln{mE / g^) + 7.0)/47r + 0{g^T). On the contrary, the vector current 
is still conserved, and thus the correlator related to VJf = ipjsip is a constant, and vanishes 
in infinite volume. Thus, the long distance screening masses related to the Ua(1) partners 
^41, V3* are different even at arbitrarily high temperatures and, in this sense, the axial Ua(1) 
symmetry is strictly speaking never restored. A similar conclusion of course also holds, for 
instance, for 5*,^*, as discussed above. 

9. Conclusions 

The purpose of this paper has been to discuss analytic predictions for the screening masses 
related to various quark bilinears at high temperatures. We have addressed both flavour 
singlets, for which the screening masses are determined by those measured with bosonic 
glueball operators in the dimensionally reduced theory, and flavour non-singlets, for which 
we have determined the next-to- leading-order perturbative correction, Eq. (|6.21|1 . 

The next-to-leading-order correction we find is small, but positive in sign. This is in 
contrast to what lattice measurements traditionally indicated, and leads to a non-monotonic 
behaviour for the mesonic screening masses, whereby they first "overshoot" the perturbative 
result 2ttT as the temperature is increased above Tc, and then approach it from above for 
T 3> Tc- It will be interesting to see whether this kind of a behaviour can be reproduced, 
as the lattice studies approach the infinite volume, continuum, and chiral limits. We should 
stress that our predictions hold formally also for the quenched theory, probably the most 
immediate case on which progress could be expected. 

Apart from 4d lattice simulations, another interesting direction for further study might be 
to implement the effective theory we have derived, NRQCD3, in lattice regularisation, and 
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measure the correlation lengths that way, following 49 (assuming that the non-trivial issues 
of renormalization and numerical instabilities in NRQCD3 can be dealt with). This would 
allow to account for some higher order corrections as well as non-perturbative effects, showing 
in particular whether those could change the sign or the magnitude of the correction we have 
discussed, at realistic temperatures not far above the deconfinement phase transition, where 
the coupling constant is not asymptotically small. 

The ultimate theoretical goal of these studies is to understand the dynamics of QCD at 
temperatures above a few hundred MeV, relevant for cosmology as well as for heavy ion colli- 
sion experiments. One may ask, in particular, whether the strict coupling constant expansion 
might work for some observables after all, when carried out to a sufficient depth (possibly 
including non-perturbative coefficients), in spite of the fact that the first few corrections are 
huge [HOI 1^ [3] . This would allow to avoid the use of more complicated and partly phe- 
nomenological tools (for which many alternatives are available, however ^H])- Apart from 
gluonic observables |12j-|17j. there are promising signs for this also from related fermionic 
observables, quark number susceptibilities [SHEl, but our case is more infrared sensitive 
and in this sense a stronger test. 
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